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Almost any textbook on the scattering theory has a chapter devoted to the Jost function, 
but none of them gives a practical recipe for its calculation. Thus one usually gets a feeling 
that the Jost function is a pure mathematical object, elegant and useful in formal theory, 
but impractical in computations. This is even more so in the case of multi-channel problems, 
where the Jost function is, in fact, a matrix. Recently, we proposed an exact method for 
direct calculation of the Jost function for central potentials (which may have Coulombic 
tails) and the Jost matrix for non-central short range potentials [^. This method 
works for all real or complex momenta of physical interest, including the spectral points 
corresponding to bound and Siegert states. In this work we extend it for potentials which 
couple channels with different thresholds. 

There are three different types of physical problems associated with the Schrodinger 
equation: bound, scattering, and resonance state problems. They differ in the boundary 
conditions imposed on the wave function at large distances. Alternatively, a solution can be 
specihed by the boundary conditions at the origin. For example, in the single-channel case 
the solution which vanishes as r —> 0 exactly like the Riccati-Bessel function, is called the 
regular solution. It is well known that if the potential couples N channels, the Schrodinger 
equation 
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has 2N linearly independent column-solutions 
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and only half of them are regular at the point r 


0. We can therefore combine all the 
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regular solutions (to distinguish them we use the notation 0) in a square matrix, 
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= ll^(J5.r)|l , 


which obeys the (matrix) Schrodinger equation with the boundary condition 
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where it{z) is the Riccati-Bessel function and /c„ = ^j2^n{,E — En) is the channel mo¬ 
mentum. Any physical solution must be regular at the origin and therefore it is a linear 
combination of the regular columns. In this respect the matrix ||0|| can be called the regular 
basis. 

In order to obtain this basis in the form which is most suitable for describing its 
long-range asymptotics, we replace the unknown matrix-function ||0(ii^,r)|| by two oth¬ 
ers, WEEI{^E, r)\\ and \\E^^\E,r)\\, via the following ansatz 
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with the additional constraint conditions 
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which make the derivative of the wave functions continuous even if the potential has a sharp 
cut-off. Substituting this ansatz into the Schrodinger equation we derive the following 
system of hrst order differential equations for these new unknown matrices. 
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The boundary conditions corresponding to (^ are 

^^^El^,{E,r) = 5nn' ■ (6) 

r—>U 

Explicit inclusion of the Riccati-Hankel functions into the basis (§) guarantees the correct 
asymptotic behaviour of the physical solution. 

In constructing a physical wave function g^niE, r) one requires the appropriate coefficients 
Cn in the linear combination of the regular columns. 


(p„(E,r) =^(j)r,n'{E,r)Cn' . 
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For bound and resonant states, for example, in the physical wave function only the term 
proportional to survives at large distances. This gives the homogeneous system of 
equations 
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which has a non-trivial solution if and only if its determinant is zero, i. e., 

det||^(-)(E,cx))|| =0 . (7) 

It is easily seen that this is the determinant of the dost matrix ||F(i7)|| which is dehned by 
the asymptotic behaviour of the regular basis ^ 

h^(~\knr)Fnn'{E) + h^g^\knr)Fnn'{E) . 

We can calculate this matrix by solving the differential equations from r = 0 up to a 
large radius where the potential is insignihcant. Since the potential vanishes, the right hand 
sides of Eqs-d^) disappear at large distances. Hence the derivatives become 

zero which means that asymptotically the matrix-functions r)|| become constants 

coinciding with ||F(i7)|| and ||F(i7)|| respectively. 

Similarly to the single-channel case (see Refs. [|l]-^) the above scheme can be easily imple¬ 
mented only for bound and scattering states. In the resonance domain of the complex energy 
plane the Riccati-Hankel functions and therefore the derivatives of ||jFR)(i7,r)|| 

at large distances, become inhnitely large which means that the limit ||jFR)(i7, cx))|| does 
not exist. Thus in order to calculate the dost matrix in this domain we need to perform a 
complex rotation of the independent variable r, r = a;exp(i6'), in Eqs.(^, where a; > 0 and 
0 < 6^ < 7r/2. This results in a movement of the unitary cuts at each threshold downwards 
as is shown in Fig. 1. After such a rotation vanishes in all areas above the turned 

cuts, which guarantees the existence of the limit cx))|| there. 
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Figure 1. Complex rotation of the co¬ 
ordinate by an angle 6 results in a ro¬ 
tation of the unitary cuts by 26. 


By locating the zeros of the dost-matrix determinant we can End the resonance energies 
and the total widths. The determination of the partial widths, however, is a difficult task 
for most methods used in the past. In our approach this is not a problem at all. Indeed, 
knowing the dost matrix, we can easily obtain the S'-matrix at a given complex energy as 

||S(B)|| = ||.Fl+l(B,«D)||.||.F'-l(B,«D)|r‘ 


which, generally, consists of two terms, viz. 
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One of them is the smooth background term and the other a singular resonance term. 
Thus a partial width is the limit 
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Practically, to find such a limit we calculate the S'-matrix at a point very close to the 
resonant energy and multiply it by the corresponding energy difference. 

To demonstrate the ability of the method we consider a numerical example proposed by 
Noro and Taylor in Ref. where they used the S'-wave potential 

which describes a two-channel system with the thresholds Ei = 0 and E 2 = 0.1 (atomic 
units). By expanding the resonance wave function in a series of square-integrable functions, 
Noro and Taylor located the hrst resonance generated by this potential. We reproduce their 
result and found also three bound states and several broad resonances. For comparison we 
put the parameters of the hrst resonance in Table 1. The digits shown do not change when 
the rotation angle is changed. This demonstrates the accuracy of our method. 


Table 1. 

Values are given in atomic units 


In conclusion, we propose a method based on simple differential equations of the hrst 
order, which can be easily solved numerically. Thus, the spectrum generated by any given 
potential can be thoroughly investigated. At the same time the physical wave function can 
be obtained having the correct asymptotic behaviour. 



our results 

Noro & Taylor 

Er 

4.7681968188 

4.7682 

T 

0.0014201920 

0.001420 

hi 

0.000051103 

0.000059 

r2 

0.001368733 

0.001361 
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